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Abstract. Here we consider stationary states for nonlinear Schrodinger equa- 
tions with symmetric double well potentials. These stationary states may 
bifurcate as the strength of the nonlinear term increases and we observe two 
different pictures depending on the value of the nonlinearity power: a simple 
pitch-fork bifurcation, and a couple of saddle points which unstable branches 
collapse in an inverse pitch-fork bifurcation. In this paper we show that in the 
semiclassical limit, or when the barrier between the two wells is large enough, 
the first kind of bifurcation always occurs when the nonlinearity power is less 
than a critical value J2}; in contrast, when the nonlinearity power is larger 
than such a critical value then we always observe the second scenario. The 
remarkable fact is that such a critical value is an universal constant in the 
sense that it does not depends on the shape of the double well potential. 



Spontaneous symmetry breaking phenomenon is a rather important effect that 
arises in a wide range of physical systems modeled by nonlinear equations. In 
classical physics spontaneous symmetry breaking occurs in optics, and it has been 
experimentally observed for laser beams in Kerr media and focusing nonlinearity [TJ 
[5] . Another natural setting in which spontaneous symmetry breaking phenomenon 
may arise is for Bose Einstein condensates with an effective double well formed 
by the combined effect of a parabolic-like trap and a periodical-like optical lattice 
[31 [H[5]. Also, the study of gases of pyramidal molecules, like the ammonia NH3, 
it is a topic where spontaneous symmetry breaking phenomenon actually plays a 
crucial role. In [6J [7] has been introduced a nonlinear mean field model of a gas 
of pyramidal molecules; in this model spontaneous symmetry breaking explaining 
the presence of two asymmetrical degenerate ground states, corresponding to the 
different localization of the molecules, has been predicted with a full agreement 
with experimental data pj [8] . 

The n-dimensional linear Schrodinger equation with a symmetric potential with 
double well shape has stationary states of a definite even and odd-parity. However, 
the introduction of a nonlinear term (which usually models, in quantum mechanics, 
an interacting many-particle system) may give rise to asymmetrical states related 
to a spontaneous symmetry breaking effect. The governing equations are nonlinear 
Schrodinger equations of Gross-Pitaevskii type 

(1) ihdtij = H q tP + e\iP\ 2 ^ , 

where e is the strength of the nonlinear term, /1 > is the nonlinearity power, and 
Hq is the linear Hamiltonian with a symmetric double well potential. When /1 = 1 
we have a cubic nonlinearity and the resulting equation has been largely studied [HI 
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[TOl \TT\ IT2l IT3|. [14] . Recently, for higher values of /i the resulting equation has been 
the object of an increasing interest with several interesting physical applications 

In the case of cubic nonlinearity a family on nonlinear stationary states bifurcates 
from the linear stationary state when the adimensional nonlinear parameter 77, 
associated with the strength e of the nonlinear perturbation by |9]). assumes the 
value rf given by equation (|12p . This nonlinear ground state branch consists of 
states having the same symmetry of the linear state, and tipically we observe also 
an exchange of the stability properties. The linear stationary state is stable for rj 
less than the value 77*, and for 77 larger than 77* the linear stationary state becomes 
unstable and the new asymmetrical states are stable: that is we have the usual 
picture of a pitch-fork bifurcation as in Figure [T] top-panel where the variable z 
belongcing to the interval [— 1,+1] represents the imbalance function. Imbalance 
function, defined in equation (|10[) . is related to the position of mean value of the 
stationary state: when z = the state in invariant (up to a phase term) with 
respect to the symmetry of the double well potential; in contrast, when z takes 
the end-point values z = ±1 then the state is fully localized inside on just one 
well (conventionally the right-hand side one for z — +1, and the left-hand side for 

* = -!)• 

However, we should remark that the picture of Figure [1] top-panel still holds true 
also for other values of the nonlinearity power, for instance for /i = 2 and /_t = 3, 
but for higher values of this parameter we observe a rather different picture [14j . 
In Figure [T] bottom-panel we consider the case of nonlinearity power /j = 5, in such 
a case a couple of new asymmetrical stationary states sharply appear as saddle 
points when 77 is equal to a given value 77 + « 4.41; then, for increasing values of 77, 
the two unstable solutions disappear at 77 = 7/* = 6.4 showing an inverse pitch-fork 
bifurcation. Thus, for 77 between the two values 77 + ss 4.41 and 77* = 6.4 we observe 
the co-existence of three stable stationary states; one of them corresponds to the 
linear stationary state which has the same symmetry properties of the potential, 
while the others two are localized on only one well. 

In this paper we investigate the bifurcation picture of the stable stationary states 
of equation for any positive value of the nonlinearity power jjL. In particular, 
in the semiclassical limit (or, equivalently, in the limit of large distance between 
the wells) we'll see that the simple pitch- fork bifurcation as in Figure [1] top-panel 
always occurs when the power fi is less than a critical value ^threshold, an d the 
couple of saddle points with an inverse pitch- fork bifurcation as if Figure [1] bottom- 
panel always appears when the power /j is larger than ^threshold- The remarkable 
fact is that such a critical value ^threshold is an universal critical power, in the sense 
that it does not depends on the shape of the double well potential and it does not 
depend on the spatial dimension. Such a critical vale is found to be equal to 

(2) ^threshold = (3 + V13)/2 . 

The linear Hamiltonian we consider 

(3) H = -—A + V, 

2m 

with a symmetric double well potential V(x) = V[S(x)], where S(x) is the symmet- 
ric spatial inversion with respect to a given hyperplane II of the Euclidean space 
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Figure 1 . In this figure wc plot the graph of the stationary states 
around the minimum of the energy (full lines represent stable sta- 
tionary states, broken lines represent unstable stationary states) 
as function of the nonlinearity parameter r\ for cubic nonlinear- 
ity (i.e. for fi = 1) in the top panel, and for higher nonlinearity 
power /j, = 5 in the bottom panel. The variable z represents the 
imbalance function. 

R". The potential has two nondegenerate minima at x = x±, x+ = S(x^), such 
that V(x) > V(x±), Vx S R" \ {z±}, and W(x±) = and Hess V(x±) > 0. If 
we consider the semiclassical limit of h small enough O [TB] , or equivalently the 
limit of large distance between the two wells [TTJ QT] , then it is well known that the 
discrete spectrum of H is given by a sequence of doublets. Let A± be a doublet 
of non-degenerate eigenvalues (A + < A_), for instance the lowest two eigenvalues 
of i?o, then there exists a positive constant C > 0, independent of h, such that 

inf [A - A±] > Ch , 
\e*(H )-{\ ± } 

where cr(i?o) is the spectrum of Ho- The splitting between the two eigenvalues 
ijj = fj(A_ — A+) exponentially vanishes as h goes to zero |18) . The normalized 
eigenvectors (p± associated to A± are even and odd real-valued functions with re- 
spect to the hyperplane II 



<p±[S{x)] = ±<p±(x). 
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The normalized right and left hand-side vectors 

tfR = (tp + + i/j_)/v2 and ip L = (ip + - </j_)/v2, 

usually named single-well states, are localized on only one well and their supports 
practically don't overlap in the sense that 

(4) max \(p R {x)ip L (x)\ = O (e~ c/h ) , as h -> , 

for some positive constant C. 

The time dynamics associated to the linear Hamiltonian ((3]) is well studied [19] : 
when the state tp is initially prepared on the space spanned by the two vectors 
Pr,l> then it performs a beating motion between the two wells with beating period 
T = 2irh/uj. Since the beating period T plays the role of the unit of time then 
we rescale the time r = u)t/H; furthermore, we consider also the gauge choice 
ip(x, t) -> e int/h ^{x, t), where ft = (A+ + A_)/2. Then equation (P) takes the form 

(5) iu)d r i> = [H -Q]ip + e\iJj\ 2f *iP 

where we apply the two-level approximation by restricting the wave-function tp to 
the space spanned by the two single well states <£r,l'- 

(6) ip = dRipR + a L ip L , 

where aR and a^ are unknown complex- valued functions depending on the time r. 
Since 

Hoij) = aR [fltpR - uj<pl} + cll [—u(fR + £1<pl] 

then, by substituting ((6]) in ([5]) and projecting the resulting equation onto the one- 
dimensional spaces spanned by the single- well states ifR and (pi,, it takes the form 
(hereafter ' denotes the derivative with respect to t) 

iu)a' R = -ua L + e(ip R , \ip\ 2 ^il)) 



(7) 



iua' L = -ua R + e((p Ll \ip\ 2 »<ip) 



where (•, •) denotes the scalar product in the Hilbcrt space L 2 (R"). From (|4|) and 
since ipR[S(x)] — (pi,(x), then a straightforward calculation led us to the following 
result 

{<p R , \M 2 ^) = c\a R \^a R + 0{e- c «' h ) 

and 

((p L , \i>\ 2 ^) = c\a L \^a L + 0{e- c -' h ) 

for some positive constants Cr and Cl, and where the constant c is the same in 
both equations and it is given by 

C = (ifR, \(p R \ 2fl ipR) = (lf L , |<^l| 2 'Vl) 

Thus, the two-level approximation ([7]) takes the following form up to an exponen- 
tially small error as h goes to zero 

/gs f ia'R = -a L + rj\a R \ 2lJ -aR 

\ ia' L = -a R + r]\a L \ 2fi a L 

where 

(9) r\ = ce/uj 
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is an adimcnsional parameter which only depends on the strength e of the nonlinear 
term and, by means of the constant c and of the splitting u), on the shape of the 
double well potential. 

We perform now the qualitative analysis of the two level approximation ([8]) by 
looking for the stationary states and studying their dynamical stability/instability 
properties. To this end we assume, for the sake of definiteness, 7/ > and let 

(10) a R = pe M , a L = qe 1 ^ z = p 2 -q 2 , 6 = a- (3 

where p and q are such that p 2 + q 2 = 1 . The imbalance function z takes value 
within the interval [—1,1] and its value is related to the interval of localization of 
the wave- function tf>. The phase 9 is a torus variable with values in the interval 
[0, 2ir). Then, ([8]) takes the Hamiltonian form 



(11) 

with Hamiltonian 



= d z n 

= -d e H 



H = 2\/l-z 2 cosO-tj 



(1 + z)» +1 + (1 - z y 

2P(p+ 1) 



which is a first integral of motion 

Equation pip always has, respectively, symmetrical and antisymmetrical station- 
ary solutions (#1, z\) and (#2, 22), where 6x = and 82 — it and where z\ — z-i = 0. 
Furthermore, asymmetrical stationary solutions may respectively occur for 9 = 
and 9 — n as solutions of equations f+(z) = and f-{z) = 0, where 



f±(z) = T2z/Vl-z 2 - rj[(l + zf - (1 - z)"]2-" . 

Since we have assumed 77 > then the derivative takes only negative values for 
any z € [— 1,+1] and thus f+(z) = has only the solution z = 0. On the other 
hand, equation f-{z) = might have other solutions coming from a pitch-fork 
bifurcation of the stationary solution z = as we can see in Figure [T] top-panel for 
•q larger than the value 

(12) 77* = lim 77(2) = 

where 



t?(z) = 2" +1 z/ ^Jl- z 2 {{\ + zf - {I- zY) 

is obtained by solving equation f-(z) = with respect to 77. 

In particular, in Figure [1] bottom-panel we observe that in the case of fi = 5 a 
couple of saddle points appears when 77 G (?7 + , 77*), where 

g(z, n) = (z 2 yL — zfJ, + 1)(1 + zY ■ 



For instance, we have that r/ + = y/27/2 w 3.67, for fj, = 4, and 77+ w 4.41 for [i = 5. 

We thus have a transition from the bifurcation picture as in the top panel of 
Figure [1] to the more complex bifurcation picture as in the bottom panel of Figure 
[TJ and the transition from the first picture to the second one occurs when the 
nonlincarity power /1 is equal to a threshold value ^threshold such that 4Up = at 
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z = 0, that is the two saddle points will merge with the stationary solution z 
Since 



d 2 r](z) 



dz 1 



2"(3^+ 1 -i?) 



=0 9 A« 

then the threshold value is given by ([2|). We may remark that this threshold in 
an universal value since it does not depend on the parameters of the double well 
model and on the spatial dimension. 

The qualitative behavior of the solutions of equation pip is then studied by 
means of the conservation of the energy 7i as done, for instance, by [6] for cubic 
nonlinearity (where [i = 1). In Figure [2] we plot the integral paths of the equation 
H = E n for some values of the energy E n , where /i = 5. In the top panel where 
T] = 2 < rj + w 4.41 we can only see closed curves corresponding to beating periodic 
motions between the two wells. In the middle panel, where rj + ss 4.41 < rj = 
5 < rj* = 6.4, we have three stable stationary solutions (circle points), two of them 
are localized on just one well and closed curves surrounding them correspond to 
periodic motions inside the well, without beating effect. In the bottom panel, where 
■q* = 6.4 < r) = 6.5; we have two stable stationary solutions (circle points) localized 
on just one well and we don't observe a beating motion around the stationary 
solution at z = 0. 



In conclusion, in this paper we have proved the existence on an universal critical 
nonlinearity power @ for nonlinear Schrodinger equations with double well poten- 
tial in the semiclassical limit. For nonlinearity power below this value we always 
observe a simple pitch-fork bifurcation phenomenon as the strength of the nonlin- 
ear perturbation increases, and the new asymmetrical stationary states gradually 
becomes localized on the single wells. In contrast, for nonlinearity power above 
([5]) we always observe a more complicate scenario: the appearance of a couple of 
saddle points where the asymmetrical unstable stationary solutions will merge with 
the stationary solution at z = drawing an inverse pitch-fork bifurcation as the 
strength of the nonlinear perturbation increases. The new physical relevant ef- 
fect associated with such a new scenario is the sharply appearance of asymmetrical 
stationary solutions fully localized on the single wells. 
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^threshold we only have the pictures of panels top and bottom, the 
picture of middle-panel does not occur. 
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